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Abstract: 

We consider the thermodynamics of a uniformly charged polyelectrolyte with harmonic bonds. For 
such a system there is at high temperatures an approximate scahng of global properties like the 
end-to-end distance and the interaction energy with the chain-length divided by the temperature. 
This scaling is broken at low temperatures by the ultraviolet divergence of the Coulomb potential. 
By introducing a renormalization of the strength of the nearest-neighbour interaction the scaling 
is restored, making possible an cfhcicnt blocking method for emulating very large polyelectrolytes 
using small systems. 

The high temperature behaviour is well reproduced by the analytical high-T expansions even for 
fairly low temperatures and system sizes. In addition, results from low-T expansions, where the 
coefficients have been computed numerically, are presented. These results approximate well the 
corresponding Monte Carlo results at realistic temperatures. 

A corresponding analysis of screened chains is performed. The situation here is complicated by the 
appearance of an additional parameter, the screening length. A window is found in parameter space, 
where scaling holds for the end-to-end distance. This window corresponds to situations where the 
range of the potential interpolates between the bond length and the size of the chain. This scaling 
behaviour, which is verified by Monte Carlo results, is consistent with Flory scaling. Also for the 
screened chain a blocking approach can be devised, that performs well for low temperatures, whereas 
the low-T expansion is inaccurate at realistic temperatures. 
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1 Introduction 



Thermodynamical properties of uniformly charged polyelectrolytes consisting of linear chains of 
monomers, with covalent harmonic bonding forces and Coulomb interactions (screened or un- 
screened), have been extensively studied e.g. with Monte Carlo methods (l| ^, H, ^, ||, ^ and 
variational techniques ^ . Less attention has been paid to the low-T and high-T expansions for 
these systems. To some extent the lack of interest in these two limits has been motivated by the 
conjecture that "realistic" temperatures fall outside the range of such expansions. 

In this work we study the scaling properties of polyelectrolytes starting out from systematic studies 
of the low-T and high-T expansions both for Coulomb and screened chains. The latter are based 
upon perturbative treatment of the interaction. Extensive Monte Carlo (MC) calculations are used 
to evaluate the results. 

For the Coulomb chain in dimensionless units there are only two parameters, the number of 
monomers N and the temperature T. For global quantities such as the end-to-end distance Tee 
and the interaction energy it turns out that one has approximate scaling for large resettled temper- 
atures T = T/N. This scaling is broken at low temperatures due to the short-distance divergence 
of the Coulomb energy. By appropriate renormalization of the harmonic and Coulomb terms and 
by introducing an additional corrective nearest-neighbor interaction, a blocking scheme results that 
allows for fairly accurate MC calculations of very large systems with quite modest computational 
investments Thus, e.g., when emulating a Af=2048 system at room temperature using a system 
of 50 blocked monomers, approximately a factor 70000 is gained in computing speed at the modest 
expense of 3.6% error in ree [f). At large N the rescaled temperature becomes small and hence one 
should be able to estimate the thermodynamical quantities using a low-T expansion. Such expan- 
sions are systematically developed, and when evaluating the coefficients numerically one is indeed 
able to compute r^e within 3% for large systems at room temperatures. 

For a Debye-Hiickel screened chain the situation is quite different. In addition to the parameters N 
and T one has the Debye screening parameter k, introducing an additional length scale for the chain. 
For large K-values the potential is short range and the chain is essentially Brownian, whereas for 
low K the Coulombic chain is recovered. As A'^ — > oo for fixed (k,T), scaling relations are obtained, 
consistent with Flory scaling When the range of the potential is much larger than the average 
bond length, but much smaller than the overall chain size, a high temperature expansion indicates 
a simple dependence on the combination N/k'^T^ of the deviation from Brownian behaviour, as 
measured by r^^/NT. The blocking approach works reasonably well also for the screened case, at 
least for relatively low temperatures, despite less firm theoretical foundation. The low-T expansions 
of global quantities, r^e and the energy, do not provide as accurate results as in the Coulomb case. 

Due to the different nature of the unscreened Coulomb and Debye screened chains our presentation 
is structured such, that each of the two systems is dealt with separately in a self-contained way. We 
also include in one of the appendices underlying details of the blocking approach of ref. |^ . 

The paper is organized as follows: In Section 2 the Coulomb model is presented together with scaling 
and scale breaking properties and numerical results. The corresponding discussion and results for 
the screened chain can be found in Section 3. In Section 4 a brief summary is given. Most of the 
details concerning high- and low-T expansions are given in the Appendix A and B respectively. 
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2 The Coulomb Chain 



2.1 The Model 

In terms of dimensionless quantities, the energy of a Coulomb chain is given by |^ 



where instead of the absolute monomer positions , the bond vectors 

= x.,+1 -Xi, i^\,...,N -I (2) 
are used, and where a runs over contiguous non-nil sub-chains, with 

r„ = ^ri (3) 

corresponding to the distance vector between the endpoints of the subchain. 

Due to the simple properties of the Hamiltonian under rescaling of the coordinates, the thermal 
expectation values for Eg and Ec are related by a simple virial identity 

2{Eg) - {Ec) = i{N - 1)T (4) 



2.2 Zero Temperature 

At zero temperature, the polymer is locked in a minimum energy configuration: an aligned config- 
uration with the nearest-neighbour distances bi satisfying 

To a good approximation, this is solved by 



N 

with C = exp(l + 7) w 4.8 (7 is Euler's constant w 0.5772). For Tee, Eq and Ec, this imphes 

Tee « 7V(l0giV)^/' (7) 

Ec ~ ^N{\ogNf^ 
Ec « NilogNf 
consistently with the virial identity (eq.(^) with T — 0. 
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2.3 Low Temperature Expansions 



At low temperature, one can make an expansion around the minimum energy configuration. Count- 
ing degrees of freedom, removing three for translational invariance and two for rotational symmetry, 
we have 

{E) - £;o + ^(3iV-5)r + 0(T2) (8) 

Using this in combination with the virial identity (eq. (Q)), we obtain for the partial energies, Ec 
and Eg, 



{Ec) 



(9) 
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{Eg) = -Eo+[-N--]T + 0{T^) 



The corresponding expansions for the correlations (r^ • Yj), needed e.g. for rge, can be found in 
Appendix A, eq. (Af;). These expansions contain coefficients that have to be evaluated numerically. 



2.4 High Temperature Expansions 

In the high-T limit, the results for and Ec can be expanded in power series in 1/T. For large 
N, the first few terms yield (see Appendix B) 

{Ec) « 

The corresponding expansion for the bond 
identity (eq.(0)). 

We note that the expressions are consistent with a simple scaling behaviour since for large N 
both Tee/N and Ec/N are functions of the rescaled temperature, T/N, only. The relevant small 
parameter of the high-T expansion is obviously N/T. In Appendix B the high-T expansion is 
discussed in some detail. 



4 /2iV5/2 
15 

4 A2"A^3/2 



3 V TT Tl/2 

energy Eg can be obtained from Ec using the virial 



2.5 Scaling Prom a Continuum Formulation 



The scaling properties indicated by the high-T expansion can be understood as follows. For large 
N , we can approximate the discrete polymer by a continuous chain, by rescaling i, and T, 

x. = A^y(^^), T^Nf (11) 

in effect replacing the discrete index « by a continuous one, t = i/N, < t < 1. Then in terms 
of y(r), the Boltzmann exponent E/T is approximated by E/T, with the effective continuum 
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Hamiltonian E given by 

E^"'^ = ['drh{T)'+ Cdr C , , , „ (12) 

h 2 Jo Jr+i/N \yir')-y{r)\ ^ ' 

The remaining A'^-dependence sits entirely in the short-distance cutoff, necessary to prevent a po- 
tential logarithmic divergence. Without the cutoff, the following naive scaling behaviour would be 
exact: 

{Ec)^Nf{f), {EG)^\N^f+]^Nf{f), {rD^N'gif) (13) 

Indeed, at high T the scale-breaking is not important; the singularity is washed out by the Brownian 
fluctuations. From the high-T expansions above, we infer the asymptotic behaviour of the scaling 
functions / and g: 

fif) « iy|f-V2 + 0(f-2) (14) 

g{f) « 3f +^A/^f-i/2 + 0(f-2) 

15 V TT 

At low T, on the other hand, the scaling is more seriously broken: 

fN{f) « {logNf'~^f (15) 



2.6 Blocking Degrees of Freedom 



Neglecting the residual iV-dependence in eq. (|12|), the scaling properties suggest the possibility of 
emulating a large system of N monomers by a blocked system oi K = N/Q < N effective monomers 
1^, using the following rescaled Hamiltonian for the latter: 

B,naive 2Q' ^ *''+^ 2^ ' ^ /-^ /-^ W^A ^ ' 

i—l i j>i ' ' 

where are distances between blocked monomers, each representing Q real monomers, while 
Q' = {N — 1)/{K — 1). This simple renormalization technique should work well at high T, where 
the detailed short-distance properties of the chain are smoothed out by large fluctuations, and less 
well at low T. 

However, the approach can be improved. The A^-dependence in the effective continuum Hamiltonian, 
eq. (p^, implies the discrepancy 

.1 /.T + l/K , I 

- ^^^^ = dr I , , ,1 (17) 

Jo Jt+1/n \y{r')~y{T)\ 

In the low-T regime, where the chain is more or less linear, the discrepancy is approximately given 
by 

E(^) - E(^) « logiN/K) f (18) 
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This indicates a possible low-T improvement to the blocking approach: simply add a corrective 
interaction term corresponding to eq. (nSf) to the blocked Hamiltonian: 



K 



= <Lve + E^(-^^+0 (19) 

Wir) = ' ' (20) 

r 

With this improved blocked Hamiltonian, very large chains can be emulated (using MC) with 
relatively small blocked systems with quite small errors all the way down to T = |^ . 

The correction term can also be derived directly for the discrete system by considering eq. (p^, with 
an unspecified corrective nearest-neighbour interaction W{r). At low temperatures we approximate 
the nearest-neighbour distance by a constant a. For the blocked system, the corresponding nearest- 
neighbour distance then is 6 = Qa. In this approximation, the Coulomb energy for large N is given 

by 

N 

TV ^ 1 N 

«7Ar«-Vy«-logiV (21) 
I 

while the blocked interaction energy for large K reads 

Uk « log K + KWib) (22) 



Choosing W{r) such that Un = Uk for b — Qa yields precisely eq. (20). 
Appendix B contains a discussion of the high-T properties of the blocked system. 



2.7 Numerical Comparisons 



In this section we compare the results from the high- and low-T expansions as well as scaling 
properties with MC data. All MC calculations, both for Coulomb and screened chains (see next 
section) were performed using the pivot algorithm |lO| ]. Most results are based on 10"* thermalization 
sweeps and 10^ measured configurations. Some of the MC results were taken from ref. j^]. For 
results from exploring the blocking method we refer the reader to ref. [p| . 



2.7.1 Lov^r T Results 



It is of interest to investigate what low T means in terms of physical temperatures. In the dimen- 
sionless units of eq. ( p]), ro om temperature (290K) corresponds to T^^O. 837808. We have computed 
ree/N , where fee=V' (^'D) using the low-T expansion corresponding to eq. ( [A5| ) for different chain 
sizes at T = Tr. In fig. |l|, comparisons are made between these low-T results and MC data. 
Also shown in fig. ^ is the T=0 solution (cf. eq. (0)). Not surprisingly the low-T approximation 
becomes more accurate as N increases since it corresponds to a lower T. For A^=1024 the error is 
less than 3%, which is quite impressive. Computing r^e for A^=1024 using eq. (A5) is 200 times 
faster than using MC. The computational demand for both methods grows like iV^. 
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Figure 1: r^e/N as a function of N for a Coulomb chain at T=Tr. Dashed line, full line and crosses 
represent T=0 (eq. (^), first order low-T expansion (eq. (|A5|)) and MC results, respectively. 

2.7.2 Scaling 



Next we test the scaling properties expected from the high temperature expansions (eq. ( |l0[ ) and 
Appendix B) against numerical results. In figs. || and |[ r^^/N and Eq /N are shown as functions 
of T = N/T. The predictions from the high-T expansion for fit extremely well for T > 0.2, 
whereas for Ec the validit y reg ion of this expansion is located at much larger T. Also shown are the 
finite- effects from eqs. (B5 B6), which are larger for Eq than for r^e as can be understood from 
eqs. (B2,B5,B6). The deviation from scaling in fig. ^ for T > 1 is almost entirely due to finite- A^ 
effects. 



3 Debye Screened Chains 



Next we turn to the case of the Coulomb chain being screened by the surrounding medium. We 
model this effect by using the Debye-Hiickel potential (being fully aware of its limitations). 



3.1 The Model 



The energy of the Debye screened chain takes the following form Q : 



E{v)=Eg + Ec = -Y.vI+Y.' 



(23) 



where the additional parameter n is the inverse screening length. In what follows we assume that k, 
is independent of T. 
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Figure 2: r^e/N as a function of T for a Coulomb chain from MC results for iV=20 (□), 40 (x), 
80 (o) and 160 (+). The lines show second order high-T expansion results for A^=20 (dotted), 80 
(dashed) and oo (full) from eqs. (B5|,B6). 
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Figure 3: Ec/N as a function of T for a Coulomb chain from MC results for A^=20 (□), 40 (x), 80 
(o) and 160 (+). The lines show leading order high-T results for A^=20 (dotted), 80 (dashed) and 
oo (full) from eqs. (pM 



In the screened case, a proper virial identity relating the average bond and interaction energies does 
not exist (cf. eq. (^). Still, one has the slightly less useful identity Q. 



2{Eg)-{Ec)-kIY,^ 



2,{N - l)T 



(24) 
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3.2 Zero Temperature 



Since the screened interaction is short-range, the ground-state configuration can, for large enough 
N, be considered translation-invariant, with a constant, iV-independent nearest-neighbour distance 
a. Defining rj = na, we have 

= - log(l - e-") (25) 

valid for Ni] ^ 1. This gives a implicitly as a function of k, and we have, for small k: 

a « (-logr/)^/^ 

■n 

(— log?7)i/3 

{-lognf^ (26) 



while for large k, we obtain 



a « ((77 + l)e-" 



1/3 



((r? + !)£-")'/' 



log K , , 

a « 3-^ (27) 

K 



Similarly, we have for the energy per monomer 



a" 



Eg/N = - (28) 
Ec/N = --log(l-e-'') (29) 



3.3 Low Temperature Expansions 



Also in the screened case, observables can be expanded in a low temperature series, analogous to the 
Coulomb case, although the details are slightly modified (cf. Appendix A). The expression for {E)^ 
eq. (^) , depends only on the nur nber of deg rees of freedom and is still valid, while the expressions for 
{Ec) and C^g) separately (eqs. |A12 , A13 ) are less transparent than the corresponding unscreened 
equations (M. 



3.4 High Temperature Expansions 

Assuming a T-independent k, and a large N ^ we have in the high-T limit 

(r^g) « 3(7V - 1)T + 4y^iV3/2K-2r-3/2 (30) 
^Ec) « (.(^^-Nn'^T^^I- (31) 
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while {Eq) in this limit is given by |(7V — 1)T + ^{Ec)- In this case, there are no obvious scaling 
properties. 



A Flory estimate of the free energy for the Debyc-Hiickel potential gives a radius of the chain 

^2/5 



(32) 



indicating that the chain will behave like a self- avoiding walk in the large- iV limit. Both Flory 
scaling and the high temperature expansion, eq. (30), are consistent with a more general scaling 
behaviour: 

/ (-^) (33) 



NT ' \k^T^ 

where the function f{x) behaves as 3 -I- 3.2 x^^^ for small x, while an asymptotic behaviour f{x) 
x^l^ would correspond to Flory scaling. 



Scaling window 



A reasonable assumption is that the relevant scaling region is where the range of the potential is 
small compared to the chain size, and large compared to the monomer-monomer separation, i.e. 

1 < K^NT < iV (34) 

In this window, the leading expression for (-Ec), eq. (pl|), has to be modified into 

{Ec) w 2Nk-^T-^ (35) 

with (Eq) = |(A^ - 1)T+ (Ec), while the expression for (rg^.), eq. (^, is unaffected. 



3.5 Blocking Degrees of Freedom 



Despite the fact that the linear chain approximation is of more limited validity in the screened case, 
one still can devise a blocking approach. Indeed, for low enough temperatures the approach works 
quite weU g. 



In this case the discrete derivation of the correction term is more convenient. Again, starting out 
from a discrete blocked energy 



E- 



(K) _ 



2—1 i j^i 



A'-l 



K 



following the same steps as in the unscreened case, one arrives at 

^2 



W{r) 



log (l - g-"'/'?) - log (1 - e""'^) 



(36) 



(37) 



In this way, the longitudinal behaviour will be rather accurately reproduced by the blocked system. 
However, at low T, the transverse fluctuations will dominate, as indicated by eq. (AS). These are 
probably more difhcult to reproduce than the longitudinal fluctuations. 
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3.6 Numerical Comparisons 



In this section we investigate the vaUdity regime of the high- and low-temperature expansions of r^e 
for the screened chain. In addition, we test the approximate scahng indicated by eq. (|3^), and the 
corresponding scale-breaking arising outside the scaling window given by eq. (^). 



In fig. ^ the high- and low-T expansions of the end-to-end distance, eqs. (30, All), are compared 
to MC data. Clearly, an increase in k enlarges the high-T interval and decreases the low-T interval. 
Quantities at room temperature T^, for biologically relevant values of k, are hard to estimate from 
both high- and low-T expansions. 



The approximate scaling as suggested by high temperature expansions and supported by Flory 
arguments (eq. (^)) is shown in fig. |5| to hold almost too good for the somewhat arbitrary cut 
Nk'^T e (8, 0.37V). In fig. | the deviation fr om this behaviour is shown by plotting h\{r1^/NT) as 
a function of hiNn^T for fixed values of N/{k'^T^). Small values of Nk^T then are equivalent to 
a high temperature Coulombic limit, while large values N) correspond to a low temperature 
Gaussian limit. Hence in both limits the curve approaches the Gaussian value, which is In 3. 

A noticeable feature of fig. ^ is the seemingly A^-independent rise to the left, indicating that in this 
region (r^g) depends only on N/ k^T^ and Nk^T, or, equivalently, on T = T/N and k = Nk. This is 
equivalent to stating that a naive (i.e. without the extra nearest-neighbour term) blocking approach 
should work well there. The scale-breaking seen in the right part of fig. where the curves bend 
downwards, occurs when N k^T becomes comparable to N , violating the assumed scaling condition, 
eq. (^H). For an infinite iV, this scale-breaking would be absent, and the corresponding curve would 
never bend down. 




-6 -4 -2 To 2 4 

ln(T) 



Figure 4: r^e/N as a function of T for a screened, iV=40, Coulomb chain with k = O.l(n) and 0.63 
(x). The lines show second order high-T expansion results (k=0.1 (dashed), 0.63 (solid)) and the 
corresponding first order low-T expansions. T^ is marked with an arrow. 
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ln(N/ {K^r) ) 
Figure 5: Inrl^/NT as a function of IuTV/k^T^ with Nk'^T e (8, 0.3N). 
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ln(N/c^T) 

Figure 6: ln(rgg/iVT) as a function of In(A^K^T). The four collection of curves corresponds to four 
different values of lnN/{K^T^): from top to bottom 5.30, 2.99, 0.69 and -3.91 for N = 40(n), 
80(x) and 160(0). The lines are only intended to guide the eye. 

4 Summary 

Scaling and scale-breaking properties of thermodynamic and conformational properties for simple 
models of charged polymer chains have been established. These properties are deduced from high- 
and low-T expansions by means of a perturbative treatment of the interaction for both Coulomb and 
Dcbyc-Hiickcl screened chains. In addition to the naive high-T limit with fixed k for the screened 
chains we also consider a restricted limit where T — > oo and k ^ while keeping Nk^T large and 
K^T small. 
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For the Coulomb chain, an increase in N makes the chain stiffer due to the long-range character 
of the electrostatic interaction. The same effect is obtained by lowering the temperature; indeed, 
thermodynamic quantities for long chains are well approximated by a low-temperature expansion. 
The relevant variable of the chain is the rescaled temperature T = T/N, as further supported by 
high-T expansions which for large N become equivalent to high-T expansions. 

For the screened chain the situation is somewhat different due to the appearance of a new length 
scale ~ the inverse screening length k~^. For K~"'^-values in a window between the bond length and 
the chain length, a relevant scaling variable N/(k'^T^) can be identified. The associated scaling 
behaviour is consistent with Flory scaling ||]. 

The scaling properties in terms of T = T/N provides a firm theoretical basis for the real space 
blocking scheme of ref. a detailed derivation is given. As reported in ref. Q this approach is 
very powerful for computing r^e and Ec, both for Coulomb and screened chains, with low errors. 
Nothing prohibits the extension of the blocking scheme to include titration and multi-chain systems. 

For relatively small temperatures, the low-T expansion from ref. is further developed into an 
efficient computational tool for computing conformational properties. In the Coulomb case results 
at room temperatures with low errors emerge from this scheme. For the screened case the low-T 
approach is not viable, at least for reasonable values of iV, T and n. 
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Appendix A. Low T Expansions 
Coulomb Chains 

This section has a certain unavoidable overlap with Appendix D in r ef. §. Low-T expansions 
for {Ec) and {Eq) are given in eq. (j9|). Computing low-T corrections to Tee is somewhat more 
elaborate. 

In terms of the zero temperature bond lengths defined by bi = W define a set of tensors 

In addition, we need the two matrices related to B, 

U = (l + 2B)-i (A3) 
V = P{1- By^P (A4) 

where P denotes the projection matrix onto the subspace orthogonal to b, which is a zero-mode of 
1 - B. 

In terms of these tensors, we have the quadratic expectation- values at low-T: 

(r, . r,) = 6,6, + T I U,, + 2V,, + + 3 ^ Ckim{b^U,k + b,U,k){Ui^ - Vim,)] + 0{T^) (A5) 

where the first two terms of the T coefficient are the naive contributions from the longitudinal and 
transverse fluctuations. The rest are corrections due to the rotational degeneracy of the T = 
configuration, which is also responsible for the transverse zero- modes (of 1 — B). In terms of this 
expression, (r^g) is obtained by summing independently over i and j, while (Ec) is obtained by 
summing over i = j and dividing by two. 



Screened Chains 



For a screened chain, with bi still denoting the (now n dependent) zero temperature bond lengths, 
the tensors take the following form 

B = E^(1 + '^M (A6) 
U = I 1 + ^2 ^ + 25 I (A7) 



V = P(\-B)^P (A8) 
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while two versions of the C-tensor are needed, one related only to parallel fluctuations C'^ppp) and 
the other, C^^**^ coupled to one parallel and two transverse fluctuations. 

(T3i,j,k ^ ^ 

In terms of these tensors we have, to first order in the low-T expansion, the expectation values 

\ ^ ^ klm / 

(All) 



cr i \ ^ ^ klm / 



~{TrU + 2TrV)+ ^^^ (A12) 



i i \ ^ ^ klm 



biUik 



(Tr t/ + 2Tr y) (A13) 
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Appendix B. High T Expansions 



Sums at Icirge N 

In what follows, we will encounter sums over subchains like the following for large N 

N 

5„ = $^L« = $^(7V-fc)fc« (Bl) 

CT fe = l 

where is the length (i.e. the number of bonds) of the subchain a. Obviously, there exist N — k 
distinct subchains of length k. 

In particular, we will encounter the following sums: 

^'^ 15 12 

S-xji = ^A^^/^ + CAT - B - ^7V-i/2 + . . . (B2) 

S_3^2 = DN - 4iVi/2 - C - ^N-^'"^ + ... 
where the numerical values of the constants A, B, etc. are 

A « -.02548, B « -.20789, C « -1.46035, D = C(3/2) « 2.612375 (B3) 
In addition, we need 

C(3) « 1.202057 (B4) 



High T expansions for a Coulomb chain 

For a system of N monomers, we have the following high-T expansions, which are obtained by 
considering the Coulomb interaction as a perturbation, 

c^. 3(iV-l)T+y|r-l/2^il/2 



At large A^ these simplify to 



riJ « 3(iV-l)T+^^r-V2>Jiy^ (B5) 

|(Ar-i)T+iJi: 



{Ea) « |(Ar-i)T+i^ir-V2^L-V2 
{Ec) « \IIt-'"Y.^-."' 



{rl;, « 3(Ar_i)T+ A./2y-i/27v5/2 (b6) 

15 V TT 

{Eg) « ^(JV-l)T+^y|T-i/2iv3/2 
{Ec) « ^y|T-V2iv3/2 
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The corresponding high-T expansions for a blocked system {N —^K — N/ Q) yield identical expres- 
sions to the order shown, as long as K is large. This can be understood by considering the blocked 
Hamiltonian of eq. (19h, with Eg defined to include the constant |(-/V — K)T, and Ec to include 



the corrective nearest- neighbour interaction. The latter term, given by eq. (|2C|), does not contribute 
in the high-T limit: it is needed only at low T. 



Generic Perturbation Expansion for Screened Chain 

For a Debye screened chain, one must distinguish between two distinct high-T limits - (1) the naive 
one, where T —>■ oo with a constant k, and (2) a modified limit, corresponding to the scaling window, 
eq. (p3[), where T oo and k ^ keeping k^T small and Nk^T large. The latter of course requires 
a large N. 

In both cases the analysis is based on a perturbative treatment of the interaction, which to a few 
orders yields the following results for a screened chain of N monomers. 

{rl) « 3{N-l)T+^K^Tj2LlF2iK^/L;f) (B7) 
{Eg) « liN^l)T+^^n'TY,L.F,{nVL;f) 
(Ec) « ^j^nY,Fi{K^Wr) 



where the Tfc() are erf-related functions, defined by 

Fk{x) = {2k - 1)!! exp(a;V2) / dy y-^'^ exp(-yV2) (B8) 



Conventional High T expansion for a Screened Chain 

In the T — > oo limit with a constant k and a fixed large iV, the sum in the expression for (y^g) is 
dominated by the large L terms (scaling) and can be approximated by an integral. In contrast, the 
sums in {Ec) and {Ec) arc dominated by small L (no scaling), and must be explicitly summed. We 
thus obtain 



{r'i^) « 3(7V - 1)T + 4y^iV^/2K"2T"^/2 (B9) 
{Eg) « \{N~l)T+^-{Ec) 
{Eg) « al)^|lNn^'T~'l^ 

The corresponding high-T expansions for a blocked system, based on eqs. (^, ^) with large K 
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and large Q, yield for (r^g) the same as above. However, (Eg) and (Ec) become too large, 



{Eg) 
{Ec) 



l{N-l)T+l{Ec) 



(BIO) 



-1/2 



c(i) 



Q 



3/2 



This is due to an overshooting by the corrective term, eq. (37), which gives rise to the Q^^^ term in 
the expression for {Ec)- On the other hand, upon omitting the corrective term, {Ec) becomes too 
small by a factor Q^^^^, due to the lack of scaling mentioned above. Obviously, for thermodynamic 
properties, the blocking method docs not apply in this limit. 



Modified High T expansion for Screened Coulomb 

With the restriction of staying within the scaling window, k^T ^ 1 <C Nk^T, all three sums scale, 
i.e. they are dominated by the large-L terms. While the expression for (r^g) remains unchanged, 
we obtain for (Ec) and (Ec) in this limit 

(Eg) « ^{N-l)T+{Ec) (Bll) 
(Ec) « 2Nk-^T-^ 



For a system blocked into K effective monomers, these expressions are reproduced, provided also 
the blocked system is in the scaling window, i.e. if k^TQ ^ 1; otherwise it yields to too small values 
for (Ec) and {Eg). 
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